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Abstract 

The frequencies and linewidths of spin waves in a two-dimensional periodic superlattice of mag- 
netic materials are found, using the Landau-Lifshitz-Gilbert equations. The form of the exchange 
field from a surface-torque- free boundary between magnetic materials is derived, and magnetic- 
material combinations are identified which produce gaps in the magnonic spectrum across the 
entire superlattice Brillouin zone for hexagonal and square-symmetry superlattices. 
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I. INTRODUCTION 



Advances in the control of spin- wave propagation and dynamics 1 have led to the demon- 
stration of magnonic bose condensation 2 and coupling of electronic spin currents to spin 
waves in hybrid systems 3 . Such effects, along with theoretical proposals to electrically- 
control spin- wave properties 4 , and theoretical suggestions of high-temperature operation 
with small switching energies, may provide the foundation for an information-processing 
technology based on spin waves^. Any such technology would benefit from magnetic ma- 
terials with designed spin-wave dispersion relations, group velocities, and linewidths. A 
common method of designing such features is the fabrication of a superlattice of different 
constituent materials, used to design electronic band structures in semiconductor superlat- 
tices, and photonic band structures in dielectric superlattices. 

Here we focus on the effect of a two-dimensional superlattice of magnetic materials on 
the magnonic frequencies and linewidths, obtained from a reciprocal-space solution to the 
Landau-Lifshitz-Gilbert (LLG) equation^. Infinite cylinders of one magnetic material are 
embedded in a second magnetic material in a periodic arrangement corresponding to a two- 
dimensional square lattice or hexagonal lattice. Large gaps within the spin- wave spectrum 
are obtained when the exchange constants and saturation magnetization of the two materials 
differ greatly; thus the gaps are considerably larger for cylinders of iron embedded within 
yttrium iron garnet (YIG) than within nickel. For iron embedded in YIG we demonstrate 
the existence of a gap throughout the superlattice Brillouin zone in the magnon spectrum for 
both square and hexagonal-symmetry magnonic crystals. In photonic crystals such a feature 
forms an essential element of photonic band gap materials 8,9 , and permits the control of 
spontaneous emission of emitters embedded within the photonic crystal; here similarly the 
spontaneous emission of magnons from a source such as a spin-torque nano-oscillator could 
be suppressed by embedding this spin-wave emitter in a fully-gapped magnonic crystal. 

Central to the accurate calculation of dispersion curves associated with a superlattice is 
the proper treatment of the boundaries between the two magnetic materials. For a magnonic 
superlattice the exchange field that enters into the LLG equations is discontinuous at the 
boundary, and that discontinuity strongly influences the spin wave dynamics. Two distinct 
forms for this exchange field have been described in the literature^ - — 2 , although to our 
knowledge it has not been pointed out that these two forms provide dramatically-different 
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solutions to the LLG equation. In Section [Til we present an explicit derivation of the correct 
form of the exchange field, followed by spin wave frequencies and linewidths for various 
magnetic material combinations in Section Hill In Section [IV] we show that solutions to 
the LLG equations for the incorrect form of the exchange field differ greatly from those for 
the correct form, and furthermore the incorrect solutions are incompatible with the spatial 
symmetry of the lattice. 



II. LLG FORMALISM FOR A QUASI-TWO-DIMENSIONAL MAGNONIC CRYS- 
TAL 

We consider a magnonic crystal composed of an array of infinitely long cylinders of 
ferromagnetic material A embedded in a second ferromagnetic material B in a square or 
hexagonal lattice; the structures are shown in Fig.fU and have lattice constant a and cylinder 
radius R cy i- The cylinders are aligned parallel to a static external magnetic field H = H z, 
and the magnetization of both materials is assumed to be parallel to H . The equation of 
motion for this system is the Landau-Lifshitz-Gilbert (LLG) equation^: 

d a (r) d 

-M (r, t) = 7/ioM (r, t) x H e// (r, t) + (r, t) x -M (r, t) . (1) 

Here 7 is the gyromagnetic ratio, M s (r) is the spontaneous magnetization, a (r) is the 
Gilbert damping parameter, and r is the three dimensional position vector. The effective 
magnetic field 

H e// (r, t) = H + h (r, t) + H ex (r, t) (2) 

acting on the magnetization M(r,i) consists of three terms: the external field H , the 
dynamic dipolar field h (r, £), and the exchange field H ex (r, t). 



Derivation of the Effective Electric Field 



We wish to derive the correct form of H e ff(r,t) to use in Eq. (CD) for our magnonic 
crystal. As shown by Gilbert 7 , the exchange field can be obtained by taking the functional 
derivative of the exchange energy. For a homogeneous material, the exchange energy is 13 

U ex [M (r)] = A| [ (Vmi (r)) 2 + (Vm y (r)) 2 + (Vm z (r)) 2 ] dr, (3) 
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where A is the exchange stiffness constant. This yields the following exchange field: 

1 SU ex [M(r)] 2A 



H ex (r) = 



r V 2 M (r) . 



(4) 



Ho SM (r) fi M* 

For the inhomogeneous crystal considered here, the values of the exchange constant and 
the spontaneous magnetization will differ for the two ferromagnets, so A and M s become 
spatially dependent quantities: 



(5) 



A(r)=A B + e(r) (A A -A B ), 
M s (r) = M SB + 6 (r) (M SA - M SB ) , 

where G (r) = 1 in material A and © (r) = in material B. The exchange energy for this 
inhomogeneous situation is 
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[M (r)] = J 



A(r) 



m x (r) 
M s (r) 



+ 



M s (r) 



+ 



M s (r) 



dr, (6) 



By approximating the energy with U ex we have neglected non-exchange terms that would 
give rise to a surface torque (such as terms in the energy associated with surface-induced 
magnetic anisotropy). 

The total magnetization will consist of both a time-dependent term and a time- 
independent term: M (r, t) = M s (r) z + m (r, t) . Using the linear magnon approximation we 
assume that the time-dependent magnetization is small compared to M s (r) and therefore 
we only keep terms up to first order in m (r, t) . With these assumptions, the inhomogeneous 
exchange field derived from Eq. ([6]) is 



2M (r, t) 



VoM s (r) 

2m(r,t) 



( V • A (r) V) 



f 



M s (r) 



( V • A (r) V) 



m (r, t) 
M s (r) 



z. 



(7) 



^oMf (r) 

The exchange field enters the LLG equation only as a cross product with the magnetization 
M(r,£). The second term is parallel to M(r, £) and thus will not contribute to Eq. (pQ). 
The third term of Eq. ((71), which is proportional to m (r, t) and parallel to M s (r), will only 
produce terms of second order in m (r, t) in Eq. (pQ) and can safely be dropped. Therefore, 
we can approximate 

H - (r - <)= ^( v -^l v ) M(r ' () ' <8) 

which produces a LLG equation from Eq. (CD) that is correct to first order in m (r, t). 



We now have the following equation for the effective field: 

H e// (r, t) = H z + h (r, *) + ^ (v • v) M (r, t) . (9) 

This form is a generalization of the boundary condition obtained at the interface between a 
ferromagnet and vacuum 14 , in the absence of any surface torque, and later derived for the 
boundary condition between dissimilar magnetic materials^^. It is also the form used in 
Ref . LI. 



B. Plane- Wave Solution to LLG Equation for Quasi- Two-Dimensional-Magnonic 
Crystal 

When solving for magnons of a specific frequency uj we write m (r, t) — m (r) exp (—iujt) 
and the dipolar field, h (r, t) = — (r) exp (—iujt), with ^ (r) the magnetostatic potential. 
With the form of the effective field in Eq. ([9]), the LLG equation (Eq. (pQ)) can be written 

iQm x (R) + M s (R) V • (Q (x) Vm y (R)) - m y (R) V • (Q (R) VM S (R)) 

-m y (R) - M * J R ) 9 * ( R ) + iQa ( R ) m (R) = 0, (10) 

y ^ Qy y 

inm y (R) - M s (R) V • (Q (R) Vm x (R)) + m x (R) V • (Q (R) VM S (R)) 

+m x (R) + M ^ q R) a ^ R) " ite* (R) ^ (R) = 0, (11) 

where = u)/(\i\ fioH ) and Q (R) = 2A (R) / (fi H M^ (R)). Additionally, since there is 
no z dependence in the above equations, the three dimensional position vector r has been 
replaced with the two dimensional position vector, R = (x, y). 

This system of equations can be efficiently solved with a plane- wave method^"—. We take 
advantage of the crystal's periodicity and use Bloch's theorem to write the magnetization 
and magnetostatic potential as an expansion of plane waves: 

m (R) = e* k R ^ m k (Gi) e* G * R , (12) 

i 

* (R) = e ikR J2 *k (Gi) e iG * R . (13) 

i 

Here represents a two dimensional reciprocal lattice vector of the crystal and k is a wave 
vector in the first Brillouin zone. The magnetostatic potential can be rewritten in terms of 
the magnetization by using one of Maxwell's equations: 

V-(h(R) + m(R)) = 0. (14) 
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Replacing h (R) with — V^(R), substituting in Eqs. (|T2l) and (fT3lh and solving for the 
potential yields 

^ (G) = z m xM (G) (G x + fc g ) + m yM (Gy + fcy) 

1 ^ (G + k) 2 " [ } 

Next, we need to be able to write the material properties M S (R), Q(R), and a (R) in 
reciprocal space. Since these have the same periodicity as the crystal lattice, this can be 
done with a Fourier series expansion: 



M S (R) = J2 M s(^)e iGvR , 

i 

Q(R) = J]Q(G,)e^- R ' 
a(R) = ^a(G0e <Gi - R . 



(16) 



The Fourier coefficients are obtained by an inverse Fourier transform: 
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M S (G) = -jAf ( (R)e- ,b 'VR. 



(17) 



where S is the area of the two-dimensional unit cell. Performing the integration for G = 
gives the average 

M s (G = 0) = M S J + M SB (1 - /) , (18) 
where / is the fractional space occupied by a cylinder in the unit cell. For G 7^ 0, we have 

,Jl (|G| R C yl) 



M s (G^0) = (M SA -M SB )2f- 



(19) 



l G l R C yl 

Here J\ is a Bessel function of the first kind, and R cy i is the radius of the cylinders. The 
following infinite system of equations in reciprocal space is obtained by substituting Eqs. 
O-flE) in Eqs. (pED and ^B>: 



itt ^2 ( m ^,k (Gi) Sij + a (Gi - Gj) m yM {Gj)) = 

3 

Y{m s (G* - G,-) (^±Ml^i±M mx k (G ) + 



(GyJ + k y ) 

H (Gj + kf 



+ ( M s (Gi - G t ) Q (Gj - Gj) ((k + G,-) • (k + G,) - (G< - G,) • (G* - G,))) 



m yik (Gj) 
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iQ ^2 ( m y* ( G i) % - « (Gi - Gj) m X;k (G,-)) = 
i 

y (m s (g, - g.) (^i±Ml^i±M m (g,) 



<L- + M s (Gi - G,) + ky) \ 
3 K jJ H (Gj + k) 2 



+ ( M s (Gi ~ Gj) Q (G, - Gj) ((k + G,-) • (k + Gi) - (G< - G,) • (G* - G,)))j m*,k (G,) }> . 

(21) 

We solve this by limiting the number of reciprocal lattice vectors in the sum and expressing 



it as a matrix equation: 



Sij 

a (Gi - Gj) 



a(Gi-Gi 



TDXX 



m = M s (Gi - G,) 



m x ^ (Gi) 

m Xtk (Gat) 
m v,k (Gi) 

jn y ^ (G N )_ 
(Gx,j + k x ) {G y j + /Cj, 



ij ij 
T>yx r>yy 



m x fr (Gi) 

m-x,k (Gat) 
«WGi) 

/%,k (G^) 



(22) 



W = + M s (Gi - Gj) 



H (Gj + k)' 



(23) 



#o(G,+k) 2 

+ J] M s (Q - Gi) Q (G t - Gj) [(k + G,) • (k + G,) - (G< - G,) • (G» - G,)fl24) 



B yx 



5 ij +M s (G i -G j ) 



{Cf x j + k x ) 
^(G.+k) 2 



+ ^ M s (Gj - Q) Q (G, - [(k + G,) • (k + G«) - (G, - G,-) • (G, - Q)p5) 
z 

The LLG equation is now reduced to finding the eigenvalues and eigenvectors for the above 
equation. 

III. RESULTS 



From Eq. f |22l ) we calculate the complex eigenvalues Q corresponding to the frequencies 
of magnons in the two-dimensional magnetic superlattices of Fe, Co, Ni, and YIG. The real 
part of fi n (k) is the magnon frequency of branch n for the wave vector k and the imaginary 
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part is the inverse spin wave lifetime. To focus on the dependence of these properties on 
magnetic material combinations we consider superlattices with a lattice constant a = lOnm, 
an external field fioHo = 0.1T, and a filling fraction / = 0.5. The material properties, M 5 , 
A, and a, are listed in Table [H 

Figures [2] and El show the empty-lattice band structures obtained from the LLG equation 
for homogeneous crystals of Fe, Co, Ni, and YIG. As the empty-lattice features are governed 
by the lattice symmetry and the material's spin wave velocity, these plots depend on material 
only in setting the frequency scale of the features. 

Figs. H] and [5] show the results for when Fe is combined with Co, Ni, or YIG. The change 
in band structure from the homogeneous case is more substantial when there is a greater 
difference in the spontaneous magnetization between the two materials. For example, the 
magnetics properties of Fe and Co are fairly similar, and so for a crystal composed of these 
materials, the band structure differs little from the homogeneous case, with only some small 
splittings of the spin wave dispersion curves occurring. However, when for a crystal of Fe 
and YIG, whose magnetizations differ by more than a factor of ten, the magnonic modes 
are almost completely different from the homogeneous crystal. Furthermore, the opening of 
band gaps in these structures is more easily attainable when the magnetization is larger in 
the cylinders than it is in the host. With Fe cylinders embedded in YIG in a square lattice, 
there are four gaps occurring within the lowest nine spin wave modes, while YIG cylinders 
in Fe shows only one small gap between the first and second spin wave modes. 

Figs. [6}|9]show the detailed dispersion curves and spin wave relaxation rates for a hexag- 
onal superlattice of Fe cylinders in Ni. Plotted are the lowest nine spin wave modes (Re (Q)) 
in the entire first Brilllouin zone as well as the corresponding inverse spin wave lifetimes 
(Im (fi)). Quality factors for these modes, corresponding to the ratio of the relaxation rate 
to the mode frequency, can exceed 100 for such spin waves, especially for the lowest-frequency 
modes. 



8 



IV. COMPARISON WITH ALTERNATE EFFECTIVE FIELD 

Some recent calculations of magnonic crystals dispersion curves used a different exchange 
fieldi£ii2 ^ an ^ e one derived i n Sec. [Til The alternate form, 



differs by the positioning of one factor of M s (r)~ outside the gradient operators. A com- 
parison of the band structures obtained for the two different exchange fields is shown in 
Fig. [TUl An examination of the band structure for a homogeneous material composed of 
Fe or Ni (Fig. [2j) indicates that the results for the derived exchanged field (Eq. ([8])) pro- 
duce a band structure that is appreciably different from the homogeneous case, whereas the 
band structures produced by the alternate exchange field ( Eq. f |26l) ) are very similar to the 
homogeneous crystal. 

In Fig. [ED we show the lowest spin wave mode and corresponding relaxation rate obtained 
for Fe cylinders in Ni when using Eq. f |26l ) as the exchange field. When looking at these 
contours, we would expect them to have the same symmetries as the real space lattice. 
For a square lattice, that would be symmetry under rotations of 90° and symmetry under 
reflections about either axis. All symmetries are present for the spin wave modes of Fig. 
[TT1 however, the spin wave lifetimes are lacking the reflection symmetries. Comparison with 
Figs. [6] and [7] shows that the use of Eq. [8] for the exchange field keeps the symmetries of 
the lattice preserved in the contours. A very slight asymmetry in Figs. [3 and [9] is caused by 
terminating the infinite summation of the reciprocal lattice vectors in the LLG equation (Eq. 
(1221) ). and disappears as the number of reciprocal lattice vectors is increased; the asymmetry 
in Fig. [H] does not. 

Thus the consequences of using the exchange field of Eq. f |26l) instead of the correct 
form, Eq. ([8]), include a dramatic underestimate of the splittings in the magnonic crystal 
dispersion relations, as well as flawed rotational symmetry of the spin wave relaxation rate. 

V. CONCLUSION 

Spin wave dispersion curves and relaxation rates have been calculated for hexagonal and 
square two-dimensional superlattices of magnetic cylinders embedded in another magnetic 




(26) 



material. The correct form of the exchange field at the boundary between these two magnetic 
materials has been found, and the difference from another form used in the literature has been 
shown to be significant. Full-zone magnonic gaps are obtained for superlattice materials that 
differ substantially in their saturation magnetization, such as Fe and YIG. Quality factors 
for spin waves can exceed 100, especially for the lowest-frequency spin mode. These results 
should assist in the design of magnonic crystals that can focus or redirect spin waves due to 
their effective band structure. 
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M s (A/m) 


A(pJ/m) 


a 


Fe 


1.711 • 10 6 


8.3 


0.0019 


Co 


1.401 • 10 6 


10.3 


0.011 


Ni 


0.485 • 10 6 


3.4 


0.064 


YIG 


0.14 • 10 6 


4.15 


0.0014 



Table I: Properties of the different materials considered for the magnonic crystals — —. 
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Figure 1: Physical structure of the magnonic crystals studied here. The ferromagnetic material B 
is the host for infinitely long cylinders of a different ferromagnetic material A arranged in either 
a square (left) or hexagonal (right) lattice. The lattice constant of the superlattice is a and the 
cylinder radius is i? cy i. 




Figure 2: Empty square lattice band structure obtained from the LLG equation for a homogeneous 
crystal of Fe (upper left), Co (upper right), Ni (lower left), and YIG (lower right) with lattice 
constant a = lOnm. Frequencies are in units of THz. 
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Figure 3: Empty hexagonal lattice band structure obtained from the LLG equation for a homoge- 
neous crystal of Fe (upper left), Co (upper right), Ni (lower left), and YIG (lower right) with lattice 
constant a — lOnm. Frequencies are in units of THz. 



14 



r 



m x r 



r 



m x r 



Figure 4: Magnonic band structures for a square lattice magnonic crystal with lattice spacing 
a — lOnm and filling fraction / = 0.5. On the left is Fe cylinders embedded in Co (top), Ni 
(middle), and YIG (bottom). The right is for an Fe host with Co (top), Ni (middle), and YIG 
(bottom) cylinders. Frequencies are in units of THz. 
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Figure 5: Magnonic band structures for a hexagonal lattice magnonic crystal with lattice spacing 
a — lOnm and filling fraction / = 0.5. On the left is Fe cylinders embedded in Co (top), Ni 
(middle), and YIG (bottom). The right is for an Fe host with Co (top), Ni (middle), and YIG 
(bottom) cylinders. Frequencies are in units of THz. 
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Figure 6: (Color online) The lowest nine spin wave frequencies (in THz) for a square lattice magnonic 
crystal composed of Fe cylinders in Ni with a filling fraction of / = 0.5, and a lattice constant of 
a = lOnm. These results were obtained using the exchange field in Eq. (j8]). 




Figure 7: (Color online) The spin wave relaxation rate (in units of GHz) corresponding to the lowest 
nine spin wave modes from Fig. El 
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Figure 8: (Color online) The lowest nine spin wave frequencies (in THz) for a hexagonal lattice 
magnonic crystal composed of Fe cylinders in Ni with a filling fraction of / = 0.5, and a lattice 
constant of a = lOnm. These results were obtained using the exchange field in Eq. (jSJ). 




Figure 9: (Color online) The spin wave relaxation rates (in units of GHz) corresponding to the 
lowest nine spin wave modes from Fig. [8] 
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Figure 10: Magnonic band structure of Fe cylinders embedded in Ni calculated from the LLG 
equation using exchange field in Eq. (j8]) (left) and Eq. (126]) (right). Frequencies are in units of 
THz. 




Figure 11: (Color online) The lowest spin wave mode (in THz) and the corresponding relaxation 
rate (in units of GHz) for a square lattice magnonic crystal composed of Fe cylinders in Ni with a 
filling fraction of / = 0.5, and a lattice constant of a — lOnm. These results were obtained using 
the exchange field in Eq. ([26]) . 
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